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Abstract 

In the large momentum transfer limit, generalized parton distributions can be calculated through 
a QCD factorization theorem which involves perturbatively-calculable hard kernels and light-cone 
parton distribution amplitudes of hadrons. We illustrate this through the Hq{x,S^,t) distribution 
for the pion and proton, presenting the hard kernels at leading order. As a result, experimental 
data on the generalized parton distributions in this regime can be used to determine the functional 
form of the parton distribution amplitudes which has thus far been quite challenging to obtain. 
Our result can also be used as a constraint in phenomenological GPD parametrizations. 
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It has been established in perturbative quantum chromo dynamics (pQCD) that hard ex- 
clusive processes in the asymptotic limit depend on the non-perturbative light-cone parton 
distribution amplitudes of hadrons|l|, i2, |3|, i^]. However, the functional form of these ampli- 
tudes in parton momenta Xi has been very difficult to determine either from experimental 
data or from theoretical calculations p. For instance, the asymptotic electromagnetic 
form factors depend only on an integral of the distribution amplitudes^, |^. In the pion- 
nucleus diffractive production of two jets, one can in principle learn about the shape of 
the quark distribution amplitude in the pion, but the process is not infrared safe as has 
been pointed out recently [7|. Poor knowledge on the distribution amplitude in the proton 
has been the main obstacle in deciding at what momentum transfer the asymptotic pQCD 
calculation is relevant [1, |^, • 

In this paper, we study the generalized parton distributions (GPDs) |Tll 1^ | of hardons 



in the large momentum transfer limit. The GPDs are a new class of hadron observables 
which combine the physics of electromagnetic form factors and Feynman parton distribu- 
tions, and are related to quantum phase-space distributions of the partons through Fourier 
transformation Apart from the renormalization scale, they depend on the momentum 
transfer t as in a form factor, light-cone momentum x as in a parton distribution, and the 
projection of the momentum transfer along the light-cone direction ^, also known as the 
skewness parameter. 

We report here that the GPDs in — t ^ oo limit are calculable through QCD factortization 
in which the non-perturbative physics is included in the light-cone distribution amplitudes 
of hadrons. Using this, the functional form of the distribution amplitudes can be studied 
through the GPDs' dependence on x and C,- Conversely, our result provide a constraint 
on phenomenological GPD parametrizations. The GPDs at large-t can be measured, for 
example, from deeply-virtual Compton scattering or hard exclusive meson production or 
doubly-virtual Compton scattering in the kinematic regime ^ —t ^ ^qcd i'^ which the 
factorization theorems for scattering amplitudes have been proven However, it can be 



experimentally challenging to measure the cross sections in this regime because of additional 
power suppression in t; we will not explore this issue here. 

We illustrate our main point first by considering the generalized parton distribution 
H{x, ^, t) for the pion, defined through 




where P and P' are the initial and final state pion momenta, respectively, t = {P — P')"^, 
and V indicates path-ordering for the light-cone gauge link. Introducing P = (P + P')/2 
along the z direction and the conjugation light-cone four-vector n, such that = and 
n - P = 1, the skewness parameter is the projection of the momentum transfer P' — P along 
P direction, ^ = —n ■ (P' — P)/2. The initial and final light-cone momenta of the quarks 
are then n ■ k = x + ^ and n ■ k' = x — C,, respectively. 

In the large momentum transfer limit, one can calculate the above GPD using the pQCD 
factorization formalism which has been widely applied to electromagnetic and other form 
factors The leading pQCD contribution is shown in Fig. 1, where the initial and 

final pion states are replaced by the light-cone Fock component with the minimal number of 
partons. The circled crosses in the diagrams represent the bilocal quark operator in Eq. (Q). 
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(b) 



(c) 



FIG. 1: Leading pQ CD diagrams contributing to the pion's generalized parton distribution H[x,^,t) 
at large —t. The circled crosses represent the non-local quark operator. 



The hard part responsible for the large momentum transfer contains a single gluon exchange 
just like in the electromagnetic form factor. In the first two diagrams (a) and (b), there is 
a hard gluon exchange between the two quark lines, and in the third one, there is a gluon 
coming from the gauge link. Since the transverse momenta of the quarks are expected to 
be on the order of Aqcd, we may ignore them in calculating the hard part. Thus we can 
effectively integrate out k± in the pion wave function to obtain the distribution amplitude, 
^^^'rijj{x, k±). The parton transverse momenta flowing into the hard part are now 



■f (27r)3 ' 



taken to be zero. 

The result of the above pQCD analysis is a factorization formula for the GPD at large t 
in terms of the quark distribution amplitude 



Hq{x,^,t,fi) = / dxidyi4>*{yi,fi)(f){xi,fi)THg{xi,yi,x,^,t,fi) 



(2) 



where T^q is the hard part and can be calculated as a perturbation series in a^. All quantities 
in the above equation depend on the renormalization scale fi. The /i-dependence in the hard 
part must be such that it accounts for the difference between the GPD and the distribution 
amplitude. 

The leading contribution to the hard part can be calculated straightforwardly. 



Tu{x,xi,yi] 



ATTagCp 



S{x-Xi] 



:i-o + 



Ai-Ai 



+ h.c. 



(3) 



where Cp = 4/3, h.c. stands for a term obtained by exchange and yi, and ^ and — ^, 
Ai = 2/1 + ^1^, Ai = Xi — Xi^ {x = 1 — x). Since < ?/i < 1, the first term contributes when 
X > C,', whereas the second term contributes when x > which indicates an up-anti-up 
pair contribution. The anti-quark is generated through the one-gluon exchange on the top 
of the valence wave function. The GPD for the down quark Hd{x, ^, t) can be obtained from 
that of the up quark through simple charge symmetry, Hd{x,^,t) = —Hu{—x,^,t). 

The above result can be translated into one for the moments of the GPDs Hjp\^,t) = 
j\ dxx''^^^Hq(x,^,t). In fact, the factorization formula applies for the individual moments, 

HjrKU) = Idx,dy,<j)*{yMx,)Ti^;i {xi,yi,^,t), where T^"g (xi, yi, ^, t) is simply the n-th 
moment of THq{xi,yi, x, ^, t). For the up quark in the pion, we have. 



n~2 



[1 - oAr^ + (1 + oAr^ + (1 - e) e ^^Ar 



■m-2 



m=0 



(4) 



which contains both even and odd powers of C,- For n = 1, the above reproduces the 
hard part in the QCD factorization formula for the pion form factor 0, 0, 13- For n = 2, 
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FIG. 2: The profile function of H{x, ^, t) for the pion in the asymptotic limit. 



TP{^,t) = ATiasCp/ixiyit) {{xi + yi + l) + 2{xi - yi)i + (xi + yi - It is easy to see 

that the linear-dependence in ^ does not contribute to H^\$^, t) because of the symmetry in 
the initial and final states. For the same reason, all odd powers of ^ in do not contribute 
to the GPD moments. 

It has been suggested from the dijet production 7] and 7*7 vr transition jl^ that 
the pion distribution amplitude at /i > 2 GeV is very close to the asymptotic amplitude 
^/QfT,x{l — x) fllQ,l3|- If so, we can make the prediction for the as follows. 



which is continuous at x = ^ and x = — ^. The quantity in the braces is a profile function 
and is plotted for four different ^ in Fig. 2. The function diverges at x = 1, indicating 
the breaking down of 1/t expansion. This divergence generates a slow decrease of the GPD 
moments at large n, and is present even when ^ = 0. We note that the limit x— > 1 and 
— t — >■ cx) may not be interchangeable. If we take the limit x — 1 first, Hg may not vanish in 
the subsequent — t ^ 00 limit because the pion momentum is now carried by a single quark. 

Now we turn to the proton case. The factorization formula for the GPD Hg takes a 
similar form 



where [dx] = dxidx2dx36{l — Xi — X2 — X3 ), and $3(xj) is the three-quark distribution 
amplitude jl7l |. 

In the leading order in a^, there are three classes of diagrams contributing to the hard 
part, each with a representative shown in Fig. 2. The first class consists of diagrams with 
two-gluon exchanges not attached to the non-local operators. Shown in Fig. 2.1 is one of the 
14 possible diagrams in this class. The second class has one-gluon coming from the gauge 
link in the non-local operator, and the third class with two gluons. To calculate the hard 
part, we arrange the first quark to have spin up, the second spin down, and the third spin 




(5) 



Hg{x,^,t) = [dx][dy]^{yl,y2,y3)^3{xl,X2,X3)THg{xi,yi,x,$,,t) 



(6) 
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FIG. 3: Representatives from three classes of QCD diagrams contributing to the proton GPD 
Hq{x,S,,t) in the asymptotic limit. 



up again, with Feynman momentum Xi, X2, for the incoming quarks, and yi, 1/2, and 
Us for outgoing quarks, respectively. We use Tj to denote the hard part with the non-local 
operator inserted on the line i. Then the hard part from the proton is 

TP = l(2Ti + T2 + T3 + T{ + T^) 

T| = ^(T2 + T3 + T^) (7) 

T/ is obtained from Tj by interchanging yi and 1/3. 
Our result for the hard part is 



t' i ^ Ai-A 



Ki3^^:^^M^^0^ + (1 ^ 3)' 
M-V2 



(A2 - ^?2)('72 - A2) J 

where = 2/3, Ai = + y^^, r/i = 1 - X3 - ?/2 + (z/2 - 2:3)^, and 772 = r]i{l ^ 2). The 
functions Kij are defined as 

i^n = l_ I L_ ^ __ 

xsVsxjyl X2y2xlyl X2y2xmxd,yi 



K12 



X2y2xlyi x^y^xly^ 

1-e 



a^22/2a;3l/3a;3 
(1 + 0(1-0 



a;22/2a;32/3 

i^2i = , (10) 

xiyix^y^x^iy^ 

and i^'j'j- = Kij{l ^ 2). T3 can be obtained from Ti by exchanging (xi,?/i) and (0:3, j/3), 
respectively. From the above, we can calculate the GPD moments for the nucleon in a 
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FIG. 4: t'^Hu{x,(^,t) for the proton at —t= 20 GeV^. CZ refers to the Chernyak and Zhitnitsky 
amplitude, GS the Gari and Stefanis amplitude, and AS the asymptotic amplitude. 

factorization form. If we take the first moment, we recover the pQCD prediction for the Dirac 
form factor Fi(Q2) QQIj. If we take the second moment, we find the pQCD prediction 
for gravitational form factors A^Q"^) and C{Q'^) at large The contribution to C((5^) 

is zero at this order in 1/t. This is because C{Q'^) also contribute to the helicity-flip GPD 
E{x,C,,t) which is sub-leading in the large —t limit. 

One can make a numerical calculation of H,iix, t) using various model amplitudes in 
the literature ^ |5|, ilfli]. Using the strategy of Ref. |l8[, we have computed t'^H^, shown 
in Fig. 4, for 3 different ^ at = —20 GeV^ with the asymptotic, Chernyak-Zhitnitsky, 
and Gari-Stefanis amplitudes. Although the CZ and GS amplitudes both give reasonable 
account of data on Ff for —t > 10 GeV^, the two yield very different predictions for the 
GPD. Note that the scale of Hu is strikingly large; a relatively small Dirac Fi is resulted 
from the cancellation in the integration. 

In summary, we have obtained a QCD factorization formula for the generalized par- 
ton distributions in terms of the non-perturbative light-cone distribution amplitudes and 
perturbatively-calculable hard kernels. We have calculated the hard kernels for the pion and 
the nucleon in the leading-order in a^. As a result, data on the GPDs in the large-t regime 
provides a way to constrain the functional form of the distribution amplitudes. 

Note added: After this paper was finished, we learned that the pion case has been studied 
in Ref. jl^, our result differs from that in the paper. Studying GPD in the large-t limit 
was first done in ji^] for 7*7 — vrvr. We thank Diehl for pointing this out to us. 

We are grateful for useful conversations with Stan Brodsky, Jian-Ping Ma, and Dieter 
Miiller. This work was supported by the U. S. Department of Energy via grant DE-FG02- 
93ER-40762. 
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